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1 
METHODS AND SYSTEMS FOR 
MEASUREMENT AND ESTIMATION OF 
NORMALIZED CONTRAST IN INFRARED 
THERMOGRAPHY 


I. CROSS REFERENCE TO RELATED 
APPLICATIONS 


This is a divisional application of application Ser. No. 
12/971,919 filed on Dec. 17, 2010, which claims the benefit 
of U.S. Provisional Application Ser. No. 61/293,426 filed on 
Jan. 8, 2010. 


II. ORIGIN OF THE INVENTION 


The invention described herein was made by an employee 
of the United States Government and may be manufactured 
and used by or for the Government of the United States of 
America for governmental purposes without the payment of 
any royalties thereon or therefor. 


III. FIELD OF THE INVENTION 


The invention is directed to the fields of non-destructive 
evaluation (NDE) and of processing of data acquired from 
an infrared camera. Specifically, the systems and methods 
described herein concern the performance of a series of steps 
for analyzing and processing digital data comprising a 
plurality of infra-red (IR) video images acquired by a system 
used for non-destructive evaluation. 


IV. BACKGROUND OF THE INVENTION 


The embodiments described herein relate generally to the 
field of non-destructive evaluation (NDE) using thermo- 
graphic images. Infrared (IR) flash (or pulsed) thermography 
is an example of a technique for NDE used in the inspection 
of thin materials such as laminated or bonded composites in 
the aerospace industry. IR flash thermography is used to 
detect delamination-like anomalies, although other anoma- 
lies, such as surface cracks, may be detected. 

An example of hardware equipment for an IR flash 
thermography system comprises a flash lamp (source of 
light/heat), a flash hood or housing, a flash power supply/ 
trigger unit, a flash duration controller, an IR camera for 
capturing video images, data acquisition electronics, and a 
computer. The computer may be used for controlling the 
flash trigger, for acquiring video data from the IR camera, 
for displaying data, and for post-processing of the acquired 
data. 

In one example of an NDE technique using IR flash 
thermography, a single sided or reflection technique is used 
wherein the flash lamp (heat source) and the IR camera 
(detector) are on the same side of a test object undergoing 
inspection. A plate is provided as the test object with a round 
delamination in the center. After applying heat to the top 
surface of the test object by triggering the flash lamp, the top 
surface area surrounding the anomaly cools faster than the 
top surface (footprint) area above the anomaly. The IR 
camera captures a sequence of images of the surface tem- 
perature in terms of pixel intensity and represents the 
anomaly as a hot spot e.g., an area warmer than the sur- 
rounding area or the reference region of interest (ROI). The 
hot spot is about the size and shape of the anomaly footprint. 
Relative pixel intensity, i.e., the difference in pixel intensity 
between the hot spot (measurement ROJ) and the surround- 
ing area (reference ROI), varies with the post-flash time. 
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Deeper anomalies appear in the IR video data at later times 
compared to the near surface anomalies. After the initial 
appearance of an anomaly in the IR video data, the relative 
pixel intensity continues to increase with time. The relative 
pixel intensity of the anomaly reaches a peak at a certain 
time, and then the relative pixel intensity decays until the 
temperature of the indication area and the temperature of the 
surrounding area become equal. 


V. SUMMARY 


In one aspect, disclosed is a method for converting an 
image contrast evolution of an object to a temperature 
contrast evolution, the method comprising: calculating a 
measurement region of interest temperature change AT; 
calculating a reference region of interest temperature change 
AT,» calculating a reflection temperature change AT,,,, 
calculating the image contrast evolution C,,/; and convert- 
ing, using a processor, the image contrast evolution to the 
temperature contrast evolution. 

In another aspect, disclosed is a method for converting a 
temperature contrast evolution of an object to an image 
contrast evolution, the method comprising: calculating a 
measurement region of interest temperature change AT; 
calculating a reference region of interest temperature change 
AT,» calculating a reflection temperature change AT, 
calculating the temperature contrast evolution C’; and con- 
verting the temperature contrast evolution to the image 
contrast evolution. 

In another aspect, disclosed is an apparatus for converting 
an image contrast evolution of an object to a temperature 
contrast evolution, the apparatus comprising one or more 
processors and one or more memory units coupled to the 
processors. The apparatus is configured and arranged to 
calculate a measurement region of interest temperature 
change AT, to calculate a reference region of interest tem- 
perature change AT,,,, to calculate a reflection temperature 
change AT,,.q, to calculate the image contrast evolution Cy’, 
and to convert the image contrast evolution to the tempera- 
ture contrast evolution. 

In another aspect, disclosed is an apparatus for converting 
a temperature contrast evolution of an object to an image 
contrast evolution, the apparatus comprising one or more 
processors and one or more memory units coupled to the 
processors. The apparatus is configured to: calculate a 
measurement region of interest temperature change AT; 
calculate a reference region of interest temperature change 
AT,» calculate a reflection temperature change AT,,4, cal- 
culate the temperature contrast evolution C’; and convert the 
temperature contrast evolution to the image contrast evolu- 
tion. 

In another aspect, disclosed is a method for assessing 
emissivity of an object being inspected by an infrared flash 
thermography system, the method comprising the steps of: 
measuring a pre-flash temperature at a measurement region 
of interest W°; calculating a camera constant C,,,,,; and 
calculating the emissivity of the object according to the 
equation of: 


In yet another aspect, disclosed is a method for calculating 
an afterglow heat flux evolution the method comprising the 
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steps of measuring a reflection temperature T,,z, and calcu- 
lating the afterglow heat flux evolution S,,,. gas, according to 
the equation of: 


4 
Sposelash OF, vel —f, Lge) 


wherein De is the reflection temperature at time of flash 
and o is the Stefan-Boltzmann Constant. 


VI. BRIEF DESCRIPTION OF THE DRAWINGS 


FIG. 1 is a schematic diagram of an exemplary single- 
sided flash thermography system in accordance with some 
embodiments described herein. 

FIGS. 2(a) and (6) are a schematic diagram of a plate 
depicting a gapping delamination with measurement and 
reference regions of interest (ROIs) and a single infrared 
(IR) video image of a test object with back drilled flat 
bottom holes, respectively, in accordance with some 
embodiments described herein. 

FIG. 3 is an infrared image of a plate with back drilled flat 
bottom slots along with foils and tape (post-flash time=0.6 
sec), in accordance with some embodiments described 
herein. 

FIG. 4 is a graph of the value of pixel intensity as a 
function of position along a horizontal cross line in a frame 
taken at 0.6 seconds after the flash, in accordance with some 
embodiments described herein. 

FIG. 5 is a graph of the pixel intensity evolution as a 
function of time on the foil ROI (lower evolution from the 
left foil ROI and upper evolution from the right foil ROI) 
with the peak time of 0.6 seconds indicated by the dashed 
vertical line, in accordance with some embodiments 
described herein. 

FIG. 6 is a graph of the pixel intensity evolution as a 
function of time for the tape ROI, in accordance with some 
embodiments described herein. 

FIG. 7 is a graph of the ROI intensity evolutions for the 
measurement (upper evolution) ROI and the reference 
(ower evolution) ROI, in accordance with some embodi- 
ments described herein. 

FIG. 8 is a graph of simple contrast evolution at the 
measurement ROI with the peak time of 0.5 seconds indi- 
cated by the dashed vertical line, in accordance with some 
embodiments described herein. 

FIG. 9 is a graph of computed evolution of the reflection 
temperature, in accordance with some embodiments 
described herein. 

FIG. 10 is a graph of estimated temperature evolutions for 
the measurement ROI (middle curve), the reference ROI 
(bottom curve), and the reflection or the background ROI 
(top curve), in accordance with some embodiments 
described herein. 

FIG. 11(a) is a graph of normalized image contrast 
computed using the pixel intensity (upper curve) and one 
fourth power of the pixel intensity, in accordance with some 
embodiments described herein. 

FIG. 11(4) is a graph of measured normalized image 
contrast (lower curve) and estimated normalized tempera- 
ture contrast (upper curve), in accordance with embodiments 
described herein. 

FIG. 12 is an exemplary screen shot displaying normal- 
ized temperature contrast evolution (upper curve) and nor- 
malized image contrast evolution (lower curve) computed 
using software that implements simulation and data analysis 
methods described herein in accordance with some embodi- 
ments. 


10 


15 


20 


25 


30 


35 


40 


45 


50 


55 


60 


65 


4 


FIG. 13 is a graph of the pixel intensity to temperature 
contrast ratio, in accordance with some embodiments 
described herein. 

FIG. 14 is a graph of simulated temperature evolutions 
with an assumed refection temperature evolution, in accor- 
dance with some embodiments described herein. 

FIG. 15 is a graph of estimated image contrast and of 
simulated temperature contrast, in accordance with some 
embodiments described herein. 

FIG. 16 is a graph of estimated afterglow flux as a 
function of time charted from the first post-flash frame, in 
accordance with some embodiments described herein. 

FIG. 17 is a graph of flash flux evolution as a function of 
time until the first post-flash frame, in accordance with some 
embodiments described herein. 

FIG. 18 is a flow diagram illustrating methods for con- 
verting an image contrast evolution of an object to a tem- 
perature contrast evolution, for modeling afterglow flux, and 
for assessing the test object emissivity, in accordance with 
some embodiments described herein. 

FIG. 19 is a flow diagram illustrating, a method for 
converting a temperature contrast evolution of an object to 
an image contrast evolution, in accordance with some 
embodiments. 

Because the embodiments described herein are subject to 
various modifications and alternative forms, it should be 
understood that the drawings and detailed description are not 
intended to limit the invention to the particular embodiments 
described herein. Instead, the following description is 
intended to cover all modifications, alternatives, and equiva- 
lents falling within the scope of the invention as defined by 
the appended claims. 


VII. DETAILED DESCRIPTION 


The objects in the drawings are not necessarily to scale 
and certain features may be shown exaggerated in scale or 
in somewhat generalized or schematic form in the interest of 
clarity and conciseness. In the description which follows 
like parts may be identified throughout the specification and 
drawings with the same reference numerals. The foregoing 
brief description of the drawings is provided for a more 
complete understanding thereof. It should be understood, 
however, that the embodiments described herein are not 
limited to the precise arrangements and configurations 
shown. Although the design and use of various embodiments 
are discussed in detail below, it should be appreciated that 
the embodiments described represent concepts that may be 
embodied in a wide variety of contexts. The specific aspects 
and embodiments described herein are merely illustrative, 
for it would be impossible or impractical to include all of the 
possible embodiments and contexts in this description of the 
invention. Upon reading this description, alternative 
embodiments will be apparent to persons of ordinary skill in 
the art. 

During the processing of digital data acquired from an 
infrared camera, normalized image (or pixel intensity) con- 
trast and normalized temperature contrast differ in value for 
objects having a value of emissivity other than one. There- 
fore, for more accurate processing of data obtained from 
objects having an emissivity value other than one, the two 
quantities should not be treated as the same. In the co- 
pending U.S. patent application Ser. No. 12/900,644, filed 
Oct. 8, 2010 and entitled “Methods and Systems for Char- 
acterization of an Anomaly Using Infrared Flash Thermog- 
raphy”, which is hereby incorporated in its entirety by 
reference, Koshti distinguishes and defines the two quanti- 
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ties of normalized pixel intensity contrast and normalized 
temperature contrast. Embodiments described herein estab- 
lish formulas to estimate the normalized temperature con- 
trast from the normalized pixel intensity or image contrast 
using the methods and systems described herein. The 
method includes steps for the measurement of the reflection 
temperature evolution by comparing the simulation tem- 
perature contrast with the measured pixel intensity contrast. 
The method also comprises the act of estimating the incident 
heat flux. Ideally, the simulation should model the com- 
pound heat source flux evolution which also includes the 
thermal afterglow. The effect of the reflection temperature in 
the pixel intensity should also be accounted for to seek a 
more refined estimation of the temperature contrast profile 
from the pixel intensity evolution data. 

The method of measuring reflection temperature evolu- 
tion is described herein. First, the reflection temperature 
evolution is established based on infrared (IR) flash ther- 
mography (or simply any transient thermography) data 
acquisition. This acquisition of data requires a test object, a 
high emissivity tape with known emissivity, and a diffused 
but highly reflective metal foil with known reflectivity. The 
method also records the steady state pre-flash temperature of 
the object using a thermocouple (or other contact sensor) or 
an accurate radiometer. The IR datacube is recorded after 
performing the technique of IR flash thermography. 

The method further comprises data acquisition from four 
regions of interests (ROI). One region is for the measure- 
ment area. The second ROI is for the reference area. The 
third ROI is for the measurement of the reflection tempera- 
ture. The fourth ROI is for the measurement of the pre-flash 
temperature of the high emissivity tape. 

Using formulas given herein, the method estimates the 
reflection temperature evolution, Then, the method com- 
putes the temperature contrast from the IR data. The emis- 
sivity factor is defined to relate the temperature contrast to 
the image contrast. 

The method uses the reflection temperature evolution to 
model the afterglow flux of the flash source. Using the 
estimated compound source evolution in the simulation 
software, the method can estimate the temperature contrast 
evolutions and then estimate the image contrast profiles on 
simulated voids. 

The invention provides an emissivity estimation tech- 
nique using the same IR camera. The technique estimates the 
emissivity of an object for a desired thermal wavelength. It 
is shown that using the foil-tape (or tape-foil) technique 
during the IR shot, the transient reflection temperature or the 
reflection temperature evolution can be recorded. If the IR 
camera is programmed to use the real-time reflection tem- 
perature using the formulas described herein, it can provide 
real-time object surface temperature. The IR camera can be 
programmed to estimate the object emissivity in real-time 
using the formulas derived herein for the foil-tape technique. 

Referring now to the drawings, FIG. 1 is a schematic 
diagram of an exemplary single-sided flash thermogaphy 
system in accordance with some embodiments described 
herein. As shown in FIG. 1, the equipment for an infrared 
(IR) flash thermography (IRFT) system 105 in accordance 
with at least one of the exemplary embodiments described 
herein comprises a flash lamp (source of light/heat) 110, a 
flash-hood 115, at least one flash power supply/trigger unit. 
120, at least one flash duration controller 125, an IR camera 
130 for capturing video images, camera data acquisition 
electronics 135, and a computer 140. The computer 140 is 
used for controlling the flash trigger unit 120 for powering 
the flash lamp 110, data acquisition of the camera video data, 
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data display, and post processing of the acquired data. The 
flash-hood 115 may be made from sheet metal and have 
highly reflective (e=0.05) inside wall surfaces. The flash- 
hood may form a box-like structure, although other struc- 
tures that form an interior volume, such as a cylinder, 
hemi-sphere, etc. may also be used. However the flash-hood 
is structurally configured, one portion of the hood has a large 
opening so as to enable the hood to be positioned over a test 
object 145. The side opposite to the opening has a hole in the 
center to provide a window for the lens of the IR camera 130 
that is positioned outside of the flash-hood 115. The IR 
camera 130 is focused at a surface of the test object 145 
positioned within the hood. At least one flash lamp 110 (two 
are shown in FIG. 1) is located within the inner wall of the 
flash-hood 115 in proximity to the IR camera 130. The flash 
lamp(s) 110 is covered with a glass shield which is highly 
transmissive for light but highly emissive for thermal radia- 
tion. The flash lamp(s) 110 when powered directs illumina- 
tion towards the test object 145 without directly shining light 
upon or into the IR camera 130. The flash-hood 115 func- 
tions as a housing that contains most of the intense flash 
from the flash lamp 110. 

If the size of the test object 145 can be accommodated 
inside the hood 115, then the test object 145 is located at the 
hood opening or slightly inside the hood 115. Otherwise, the 
test object 145 is located slightly outside of the hood 115 
opening. A short duration (e.g. 3 millisecond), intense (12 
kJ) flash is triggered by the computer 140. Data acquisition 
is initiated a few seconds before the flash, and it continues 
until a prescribed duration of time has expired. The camera 
130 provides a sequence of IR images (or frames) called the 
data-cube (or digital video) of the test object surface taken 
at the chosen frame rate (e.g. 60 Hz or 60 frames per see). 
As described in further detail herein, the intensity (numerical 
value) of each pixel of the image is a function of the surface 
temperature of the corresponding area observed of the test 
object 145 at the time of the image frame. The flash causes 
the surface of the test object 145 to warm up slightly, and the 
heat then starts to dissipate rapidly. The surface cools 
through thermal radiation, convection, and conduction. Heat 
conduction within the test object 145 is considered to be the 
dominant heat transfer mode until the temperature gradients 
within the test object become small. At later times, the heat 
conduction is of the order of the combined effect of heat 
convection and radiation. The IR data acquisition and data 
analysis utilizes the thermal data in the short duration 
immediately after the flash where the thermal dissipation is 
dominated by the heat conduction within the test object 145. 

Heat exchange across the boundaries due to convection 
can be assumed to be zero if the Biot number (Nz,=hL/k) is 
less than 0.1. For example, a half centimeter thick graphite/ 
epoxy (k=0.64 W/mK) plate, using h=10 W/m°K, has a Biot 
number of 0.078. Therefore, heat conduction is the dominant 
mode of heat transfer in this graphite/epoxy plate example. 
Thinner test objects tend to equalize the temperature within 
the test object very quickly and have a relatively longer 
cooling time by heat loss to environment. 

FIG. 2(a) is a schematic diagram of a plate serving as the 
test object 145. FIG. 2(a) depicts a gapping delamination or 
anomaly 220 with measurement and reference regions of 
interest (ROIs) in accordance with embodiments described 
herein. FIG. 2(6) shows a single infrared (IR) video image 
of a test object with a flat bottom hole in accordance with at 
least one embodiment. Referring now to FIGS. 1 and 2(a), 
the test object 145 is made of a thermally isotropic material 
with constant thickness 205 that fits inside the hood 115. The 
test object 145 is supported at the corners on insulating 
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standoffs 150 and the hood is oriented vertically thereabove. 
Assuming that the flash intensity is uniform over the test 
object top surface 210, the heat conduction will be in a 
direction normal to the test object top surface 210 in most of 
the acreage area (i.e., area away from edges of the test object 
and flash boundary). The heat is conducted uniformly from 
the top surface of the test object 210 to the bottom surface 
of the test object 215. Normal heat conduction will be 
obstructed by an anomaly such as a small round gapping 
delamination 220 at the center of the test object, as depicted 
in FIG. 2(a). The delamination 220 is at some depth (d) 221 
below the top surface of the test object 210. The delamina- 
tion may also be characterized as having sonic diameter (D) 
or width 222 and a gap 223. The volume bounded by the 
anomaly 220 on one side and the top surface 210 on the other 
side is called the heat trapping volume 225. 

The top surface area 230 surrounding the anomaly 220 
cools faster than the top surface (footprint) area 235 above 
the anomaly 220. The IR camera 130 captures the surface 
temperature image in terms of the pixel intensity and shows 
the anomaly 220 as a hot spot (e.g. an area warmer than the 
surrounding area) which is about the size and shape of the 
anomaly footprint. The relative pixel intensity of the hot spot 
changes with the time. Deeper anomalies appear at later 
times in the IR video compared to the near surface anoma- 
lies. After the appearance of an anomaly in the IR video, its 
relative pixel intensity continues to increase with time. The 
relative pixel intensity associated with the anomaly reaches 
a peak at a certain time, and then the relative intensity decays 
until the indication area temperature and the surrounding 
area temperature become equal. 

FIG. 2(6) shows one of the frames in a sample IR 
data-cube from a reinforced carbon-carbon (RCC) part with 
round flat bottom holes (FBHs) machined into the back 
surface to simulate gapping delaminations of differing 
depths. IR indications 240 from different flat bottom holes 
are clearly visible in the IR image. A measurement ROI 245 
is shown as a square in the middle of one of the anomalies. 
A reference ROI 250 is shown as a square outside the 
anomaly area. As the test concludes, the test object continues 
to cool down to ambient temperature through heat convec- 
tion and radiation. 

In some embodiments, such as for a test object that is a flat 
plate with a large diameter, the top surface may be repre- 
sented as being at a temperature T, and similarly the bottom 
surface is at a temperature T,. A heat conduction rate for an 
isotropic slab may accordingly be given by equations (1) and 


(2), 


2 fpsthe 
~ (LTkA) 


(1) 


qk 


L (2) 
R= kA’ 
where q, is the heat conduction rate (cal-sec™' or BTU-hr“'), 
k is the thermal conductivity (cal-sec"'-° C.-'-cm™! or 
BTU-hr“'-° F.-!-ft"'), A is the cross sectional area of the 
plate (cm? or ft”), L is the thickness of the plate (cm or ft), 
and R, is the conductive resistance (sec-cal7'-° C. or 
hr-BTU"'!-° F.). Equations (1) and (2) may be used for 
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representing heat conduction across a delamination gap L 
with one surface at temperature T, and the other at tempera- 
ture T,. 

In some embodiments, such as for a flat plate with a large 
diameter, the heat convection rate from a surface is given by 
equations (3) and (4), 


Tsurg — Tam (3) 
= ———— and 
caer (0) | ae 
1 (4) 
al Sr 


where T,,,,¢ is the top surface temperature, T,,,,, is the 
ambient temperature of the air contacting the surface, q,, is 
the heat convection rate (cal-sec™! or BTU-hr7’), h,,,,,, is the 
heat convection film coefficient (cal-sec™'-° C.-'-cm™?, 
BTU-hr! F7!-ft-? or W-K7'-m7?), A is the surface area of 
the plate(cm”, m? or ft*), and R,, is the convective resistance 
(sec-cal“'-° C. or hr-BTU™!-° F.). For natural (unforced) 
convection, a typical value for the heat convection film 
coefficient is 4 W-K~'-m~?. Using equations (3) and (4), the 
heat convection rate across each side of a delamination air 
gap with one surface at temperature T, and other surface at 
temperature T, (T,>T,) is given by equation (5a), 


a T, — Tamb = Tomy — Tr 
(L/FeonvA) — (L/ BeonvA) 


(Sa) 


Gh 


Further, the ambient temperature of the air for thin gaps is 
given by equation (5b), 


T, +T) 
2 


(5b) 


Tamb = 


When the air gap becomes smaller than a threshold value 
represented as the Biot number, the amount of heat exchange 
by convection is replaced by heat conduction. The heat 
exchange across the boundaries due to the convection can be 
assumed to be zero if the Biot number is less than or equal 
to 0.1 (Nz,=hL/ks0.1). This relationship gives an expression 
for the threshold air gap thickness as L,,,,=0.1 k/h. Taking the 
air conductivity to be 0.026 W/mK and assuming the con- 
vection coefficient to be 5 W/m?K, then a threshold gap 
thickness of 0.020 inch (0.5 mm) is calculated, beyond 
which heat convection is applicable. Heat conduction is 
present if the gap thickness is less than the threshold gap 
thickness. The narrowing of the air gap results in increased 
heat conduction across the air gap, which manifests as a 
decrease in the peak contrast value of the void. Loss of the 
peak contrast value by 20% may not occur until the gap 
narrows significantly compared to the threshold value of the 
gap thickness. 

In some embodiments, such as for a flat plate with a large 
diameter, the radiative heat transfer rate is given by equa- 
tions (6), (7), and (8), 


Tsu — Tren (6) 


 CThA) 
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1, =€O( East Trop) Curg +E)» ad ) 


1 (8) 


where T,7 is the apparent temperature experienced by the 
object surface due to the camera-side background tempera- 
ture in units of Kelvin or Rankine, T,,,,,, is the temperature 
of the air surrounding the object in Kelvin or Rankine, o is 
the Stefan-Boltzmann. Constant 5.670400x107* W-m-?-K~*, 
€ is the emissivity of the test object surface (assuming same 
emissivity on top (front) and bottom (rear) surfaces), h,.is the 
heat radiation transfer coefficient (cal-sec7!-° C.~'-cm™? or 
BTU-hr"!-° F.-!-ft-?), and R,. is the radiative resistance 
(sec-cal7'-° C. or hr-BTU7'-° F.). At a temperature of 300 K 
for a 10 K temperature difference between the object surface 
temperature and the corresponding reflection (background) 
temperature, the radiation transfer coefficient is less than 6 
W-K7l-m™?. 

In some embodiments, such as for a large but thin slab of 
isotropic material with uniform thickness that is exposed to 
a heat pulse on the front surface, the resulting one dimen- 
sional transient heat transfer is described by equations (9) 
and (10a), 


aT aT 4 (9) 
We ae 
T(=0)=T yp, (10a) 


where «. is the test object (slab) diffusivity (cm?/s), and z is 
the distance along the thickness direction. The total heat 
transfer rate at the slab surface is given by equation (10b), 


GUAT d- (10b) 


Neglecting the air emission, the heat flux boundary condi- 
tions are given by, equations (11) and (12), 


-{k- 


aT 
ad 


(11) 
} = 9-H front T front — Tamb—front) — TAT pont — Trept— front) 
front 


aT 
ad 


Pm (12) 


refl-rear)» 


{k- 


where T(t=0)=T,,, is the initial steady state temperature of 
the test object before the flash and h,,,, and h,.,, are the 
convection coeflicients at the front and rear surfaces, respec- 
tively. 

In practice, an IR camera measures pixel intensity. Radio- 
metric cameras can measure the surface temperature. Radio- 
metric cameras are not known to be set-up for compensating 
the transient surrounding temperature. Radiometric cameras 
may have input for a constant reflection (background) tem- 
perature. They may use a reflective foil to assess the stable 
reflection temperature. For the exemplary embodiments 
described herein, either a radiometric IR camera (i.e. mea- 
sure apparent temperature) or an IR camera that measures 
the irradiance (i.e. measure pixel intensity) may be used 
because the exemplary embodiments use a method described 
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herein to capture the transient reflection temperature. The 
method assumes that camera pixel intensity response is 
proportional to the camera irradiance in the desired tem- 
perature range. 

Considering the definition of the normalized contrast 
based on the pixel intensity, its relationship with the tem- 
perature of the test object will be derived. The pixel refers 
to the camera image picture element. The pixel size refers to 
the instantaneous field of view (IFOV), which is the corre- 
sponding area of the object surface imaged in the pixel. A 
selected pixel grid area, usually in the form of a rectangle, 
is called the image region of interest (ROI, or when plural 
ROIs). The size of the ROI then refers to the area of the 
object surface imaged in the ROI. The corresponding area of 
the test object is referred to as the object ROI. The camera 
measures total radiation incident on its detector array ele- 
ment and displays it as a pixel in the image. The detector 
element response or the pixel intensity registered by the 
camera is governed by the following equation (13a), 


W= Wout WeeptWV oom 


obj 


(13a) 


where W is the average pixel intensity due to the heat 
irradiance (measured in units of W/m?, cal-sec™'-cm™, or 
BTU-hr"'-ft-?) of the corresponding camera detector ele- 
ments from the object ROI measured in gray scale bit value 
(positive integers). The right hand terms in equation (13a) 
are given by equations (14), (15), and (16), 


W, 


‘obj EUW" 


obj? (14) 


Wreg(1-eeW'yoq, and (15) 


Wom= UAT) coms (16) 


where W,,,, is the contribution to the pixel intensity due to 
the heat emission from the object ROI, W,, is the contri- 
bution to the pixel intensity due to the heat reflection from 
the object ROI, and W,,,,,, is the contribution to the pixel 
intensity due to the heat emission from the air between the 
object ROI and camera. If the method assumes a perfect 
focus and that the ROI has a uniform temperature at any 
time, then equations (17), (18), and (19) will apply, 


W'o0j=Coam0T', 


cam 


(17) 


W'rep=CoamOE ref, and 


(18) 


_— 
W atm 


Cc. 


cam 


of, 


cat 9 (19) 
where W',,, is the uncompensated contribution to the pixel 
intensity due to the heat emission from the test object ROI, 
W',q is the uncompensated contribution to the pixel inten- 
sity due to the heat reflection from the test object ROI, Wy. 
is the uncompensated contribution to the pixel intensity due 
to the heat emission from the air between the camera and 
object, T is the surface temperature in Kelvin or Rankine, 
T,.q 18 the apparent temperature experienced by the object 
surface due to the camera-side background temperature in 
Kelvin or Rankine, T,,,,, is the temperature of the air between 
the camera and the object in Kelvin or Rankine, C,,,,,, 1s the 
camera constant, € is the thermal emissivity of the test object 
surface, and t is the thermal radiation (at camera detection 
wavelength) transmissivity of the air between the test object 
and the camera. The thermal emissivity is dependent on the 
thermal wavelength for a non-gray body. The wavelength 
spectrum is affected by the surface temperature. Assuming 
that the emissivity is high and the test object is a gray body, 
the emissivity is not dependent on the wavelength. However, 
the IR camera measures radiation in a selected band of the 
wavelength (e.g. 3-5 micron). Equations (17), (18), and (19) 
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apply to the total radiation integrated over all wavelengths. 
Thus, in reality, in equations (17), (18), and (19), the 
exponent (or power) of the temperature variable T may not 
be equal to four. The exponent of the temperature variable 
may be between 2 to 17. The following analysis is applicable 
for these values of the exponent. 

A calibrated camera provides a linear response with the 
thermal radiation incident on the detector array element. 
This relationship implies that, in a hypothetical case of 
absolute zero temperature, if the camera calibration is 
extrapolated, all detector elements shall provide zero pixel 
intensity. The method assumes that the camera constant is 
the same within a range of 280 K to 350 K. The pixel 
intensity or the pixel grayscale value also depends upon the 
pixel bit resolution (typically 14 bit). The pixel intensity has 
an upper saturation limit dictated by the pixel bit resolution. 
Atypical InSb mid-wave camera used in flash thermography 
operates at a wavelength of about 3 to about 5 micron. The 
transmissivity of air is a function of the distance, humidity, 
air composition, and wavelength. A pronounced dip in the 
transmissivity occurs at a wavelength of 4.2 micron due to 
CO, absorption. Due to a short distance (~1 ft) between the 
camera and the test object, the contribution of air emission 
(emissivity~0) is very small and is neglected herein. Cor- 
respondingly, the air transmission (transmissivity~1) is 
almost 100%. The microbolometer long-wave camera works 
in the 9 to 12 micron wavelength and has better air trans- 
missivity even for longer distances of several feet (e.g. 30 
ft). Therefore, equation (13a) can be simplified to equation 
(13b), 


WaW. 


objt Weg (13b) 


The normalized image contrast is defined based on the 
pixel intensity as expressed in equations (20), (21), and (22), 


a. AW —AWyer (20) 
Wo= Aur laa 

AW + AWyer 
AW=W-W, and (21) 
AW Wrap Wreps (22) 


where C,,’ is the normalized IR image (pixel intensity) 
contrast, AW is the change in the pixel intensity of the 
measurement ROI after the flash, W is the pixel intensity at 
the measurement ROI at the post-flash time t, W° is the pixel 
intensity at the measurement ROI before the flash, AW,,, is 
the change in the pixel intensity of the reference ROI after 
the flash, W,,.r is the pixel intensity at the reference ROI at 
the post-flash time t, and Wee is the pixel intensity at the 
reference ROI before the flash. The normalized image 
contrast is also called the normalized irradiance or pixel 
intensity contrast. 

The normalized surface temperature contrast is defined as 
expressed in equations (23), (24), and (25), 


oe _ AT AT reg (23) 
AT +AT, ef” 

AD of=Trep-T reps and (24) 

AT=T-T°, (25) 
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where C’ is the normalized flash surface temperature con- 
trast, AT is the change in the surface temperature of the 
measurement ROI after the flash, T is the surface tempera- 
ture at the measurement ROI at the post-flash time t, T° is the 
surface temperature at the measurement ROI before the 
flash, AT, the change in the surface temperature of the 
reference RO] after the flash, T,,,is the surface temperature 
at the reference ROI at the post-flash time t, and Tee is the 
surface temperature at the reference ROI before the flash. 

The normalized image contrast is related to the normal- 
ized temperature contrast. As a first approximation, an 
exemplary embodiment of the methods described herein 
assumes that the air transmissivity is unity due to short 
distance (e.g. 1 ft) and the air heat radiation is negligible 
between the object and the camera. Therefore, equation 
(13a) simplifies to equation (26), 


WezeW' 


‘nt (1-€) Weg (26) 


Substituting equations (17) and (18) into equation (26) 
results in equation (27), 


WeeC camO(€T'+(1-€) Drog')- (27) 


Assuming that the contribution to the pixel intensity from 
the reflected temperature remains the same between the 
measurement ROI and the reference ROI, the temperature 
rise at the measurement (object) ROI is related to the rise in 
the pixel intensity at the corresponding image ROI by 
equations (28) and (29), 


AWC gan (TIT 4(1-€) Lei! Trep)), and (28) 
AW pe®CcamO(€(Tyep! Trop? +(1-€)(ZrepTrop? (29) 


Substituting the equations (28) and (29) in the definition of 
image contrast, equation (20), the method arrives at the 
expression in equation (30), 


o((7* - 7) - (14, - 725) ce 


Cy = 
WT (ols =) + (Thy = 1%)) +2 - 2) Ty — Ta) 


If a simulation program, such as the software program sold 
commercially known as ThermoCalc, is used to generate the 
temperature evolution at the measurement and reference 
ROIs, then equation (30) can be used to generate the image 
contrast evolution provided the method estimates the reflec- 
tion temperature evolution. In practice, the reflection tem- 
perature changes with time due to cooling of the flash lamps. 
The reflected temperature evolution can be estimated based 
on the measured data on a high reflectivity foil and a high 
emissivity tape. With the estimation of the reflected tem- 
perature evolution, the expression of equation (30) can be 
used to estimate the normalized image contrast. 

The estimated image contrast would be partially compen- 
sated for the reflection temperature. Only the camera mea- 
surement for the reflection temperature is compensated. The 
reflection temperature evolution is a measure of the heat 
source (or incident) temperature evolution. The net heat 
transfer due to the radiation is relatively small. Heat con- 
vection is accounted for in the simulation performed by the 
ThermoCale software. The simulation program provides 
choices for the shape and duration of the heat pulse. Depend- 
ing upon the user’s choice, the heat source simulation may 
or may not model the afterglow or the reflection tempera- 
ture. Thus, if the reflection temperature is not modeled in the 
simulation, then the estimation of the image contrast from 
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the simulated temperature would be considered to be par- 
tially compensated for the afterglow. 

Alternatively, a second approach is to increase the con- 
vective heat transfer coefficient by the average radiation heat 
transfer coefficient. The convective heat transfer coefficient 
is of the order of 5 W/m?K for flash thermography. 
Although, the radiative heat transfer coefficient is a function 
of the temperature, it is less than 5 W/m?K for a 10 K 
difference in the temperature range of about 300 K ambient 
temperature. The radiative heat transfer coefficient can be 
lumped together with the convective heat transfer coefficient 
to a value of 10 W/m?K in the simulation. This relationship 
may allow better prediction of the surface temperature in the 
flash thermography using the model and also better predic- 
tion of the image contrast. 

Yet another approach is to model the heat source as a high 
intensity pulse followed by a very low intensity decaying 
afterglow. 

The estimation of the temperature contrast from the IR 
measurement requires measurement of the pre-flash surface 
temperature and the measurement of pixel intensity of a 
black tape affixed to the surface of the test object. For 
instance, the method may use a high emissivity (€>0.95) 
tape (e.g., tape sold commercially by the 3M Company 
under the product name of 3M 33 or 3M Scotch Brand 235) 
positioned near the test object. The emissivity of the tape is 
assumed to be constant for both the pre-flash and the 
post-flash times. The method measures the pre-flash surface 
temperature of an ROI in the center of the black tape using 
a calibrated radiometric IR camera. The surface temperature 
needs to be close to the equilibrium or steady state. The 
method also measures the ROI intensity and assumes a high 
transmissivity value for the air. Thus the pixel intensity due 
to the tape can be expressed as given in equation (31), 


Wrape=CeanO((1-€ ape) Dre +€tapel rape’): (31) 


Due to the high emissivity of the tape, equation (31) can be 
simplified to equation (32), 


W, 


abe 4 
tape Coan DEapel; tape ) . 


(32) 


If the method measures the pre-flash tape temperature, then 


the method can estimate the camera constant (C',,,,,)as given 
in the expression of equation (33), 
We (33) 
Clam = CeamO = —P.. 
ftape Tepe 


Later a more precise expression for the camera constant will 
be derived. 

To measure the T,,, or W,,.,, the method introduces a 
diffused high reflectivity (€<0.05) reflector such as an alu- 
minum foil that is adequately wrinkled. The emissivity of 
the foil is assumed to be constant for both pre-flash and 
post-flash times. The foil is placed level with the object and 
close to the physical location of the measurement and 
reference ROIs, such as on the top surface of the test object 
surface. The method can use two foils, one on either side of 
the measurement ROI. Here, the method assumes that the 
reflection temperature is uniform over the test object and the 
foil, as given by equation (34), 


W pot C camo (1 -€yoi) T, ef toil pout )- (34) 


The contribution of the foil emission is relatively small 
compared to the foil reflection. In the exemplary method, 
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from equation (34), an expression for the reflection tem- 
perature may be written as expressed in equation (35), 


(35) 


_ wo E foil Toit 
tape| Erape pot 


tape 


Cram(1 — €foit) 


Also, the exemplary method assumes that the foil tempera- 
ture is about same as the pre-flash reference temperature due 
to the contact with the test object. The black tape is in 
contact with the test object, therefore the temperature esti- 
mations of equation (36) hold 


Djoi® Tpoit = Trey ST ape (36) 
Equation (35) then simplifies to equation (37a), 
Epi (37a) 
W pot — Wipe = ! 
TA a tape 
i Cram(1 - €foit) 


The pre-flash reflection temperature may then be expressed 
by equation (37b), 


Efe if (37b) 
Whit a wi. a 
Tren ¥ 
on Clan (1 = €foit) 


The difference between the fourth power of the reflection 
temperature before and after the flash is now given by 
equation (38), 


0 
TY = Writ — Wroit (38) 


7m 7 Clam(l = € it)” 


Knowing the camera constant, equation (38) can be evalu- 
ated as a function of the post-flash time and then substituted 
in the image contrast expression of equation (30). The 
combined steps of using the high reflectivity foil and of 
using the high emissivity tape is referred to as the “Foil-Tape 
Method” for the exemplary embodiments described herein. 


The surface temperature evolution data can be obtained 
using a simulation program such as the previously men- 
tioned ThermoCale software. The IRFT data acquisition 
using a radiometric camera that accounts for the transient 
reflection temperature can provide the surface temperature 
evolution data. Currently, this feature is not available in IR 
cameras. Here, the exemplary embodiments provide an IR 
measurement method so that measurement of the transient 
surface temperature can be performed by the IR radiometric 
cameras. The change in the pixel intensity at the foil ROI 
may be defined as expressed in equation (39), 

AW joi W foit- Wroit- (39) 
Using equations (30), (38), and (39), the estimated image 
contrast expression using the highly reflective foil is given 
by equation (40), 
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A simpler expression for the normalized image contrast will 
now be derived. From equation (13a), the change in the 
measurement RO] intensity may be derived as expressed in 
equation (41), 


AW = Ceamo[e(T - TT + 797? +1”) + (41) 


Ut ieee gr Ceres Gan Vaart pan 


Similarly, the change in the intensity for the reference ROI 
may be expressed as in equation (42), 


AWrep = Coamo [Tre — Top Treg + Top (Tap + Top) + ©) 


(1 -8\(Trep — Tog Tren + T2q)(T2q + T%))- 


The difference between the change in the measurement ROI 
intensity and the change in the reference ROI intensity may 
be expressed as in equation (43), 


AW = AWyep = Coamo[e(T = TNT +1°(T? +7) - (43) 


2 
(Treg — Trop (Treg + Trop (Teg + Trop) 


The sum of the change in the measurement ROI intensity 
and the change in the reference ROI intensity may be 
expressed as in equation (44), 


AW 4 AWrer = Ceamofe(T — TOT + TT? 47°) + (44) 
f 
2 
&(T reg — Trop (Treg + Trop )( Treg + Trey) + 


2 
21 = EN Tren — Teg) Tren + Tep\(Taan + Ten))- 


Assuming that the difference between the reference and 
measurement temperature is small (e.g., <5 K) at the time of 
maximum contrast, then the expression in equation (45) may 
be shown, 


(T+ 7 (72 +7) = (45) 


2 2 
(Treg + Trop (Treg + They) = Tren + Treg) Trent + Tren) 


Therefore, using equations (43) and (44) the expression in 
equation (46) may be shown, 


AW -AWyep [ef - T°) - eyes — Ths) (46) 


AW + AW yep [e(T — T°) + (Tre - Tey) +20 -2)Tren -— Teg) 


Using equation (20), equation (46) may be written as 
expressed in equation (47), 
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(AT -AT,.f) 
© [AT + AT ep) +20. — AT ea)’ 


(47) 


Cw 


where the change in reflection temperature is expressed in 
equation (48), 
AT, ¢-(Qef, eft): 


Therefore, equation (47) may be written as expressed in 
equation (49), 


(48) 


(AT — AT yep) (49) 


ae 


T % 
(AT + ATref) + - - 1}a7a| 


Using the definition of the temperature contrast of equation 
(23), the method arrives at the normalized image contrast 
expression of equation (50), 


Cc (50) 
om 


I 
2 - . 1a 

1 Se 
+ TAT + AT np) 


Next, equation (50) may be rearranged for the normalized 
temperature contrast expressed in equation (51), 


ze 


Thus, the normalized image (pixel intensity) contrast and the 
normalized temperature contrast can be approximately 
related to each other. An expression for the normalized pixel 
intensity contrast to normalized temperature (or pixel inten- 
sity/temperature) contrast ratio e' may be expressed as given 
in equation (52), 


1 (51) 
{= 7 1aTen 
(AT + AT ye) 


; 1 (52) 


ge = ——_________., 


I 
|; (= - 1a 


+ TAT + AT) 


and an expression for the relationship between the normal- 
ized image contrast and the normalized temperature contrast 
may be written as expressed in equation (53), 


Cy'aeT. (53) 


The pixel/temperature contrast ratio indicates that the reflec- 
tion temperature rise (caused by flash decay), the object 
temperature rise (due to the flash power during rise time), 
and the emissivity affect the conversion between the two 
contrasts. The reflection temperature is due to the glass 
shield and the flash lamps. An estimated amount of less than 
50% of the flash energy comes through the glass shield of 
the flash lamps during the flash duration. The remaining 
energy comes through the glass shield as the flash lamps 
continue to glow after the flash, providing an increase in the 
reflection temperature. Similarly, the, glass casing of the 
bulb and the glass shield get hot during the flash and 
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continue to provide thermal radiation. Some flash-hood 
designs use fans to cool the lamps after the flash to reduce 
the amount of heat from afterglow. Assuming that the 
reflection temperature is the same as the ambient tempera- 
ture immediately after the flash, then equation (53) may be 
simplified to the expression in equation (54), 


CeCy. (54) 


A shutter that covers the flash lamps immediately after the 
flash can be used to achieve a constant reflection tempera- 
ture. In practice, AT, is positive and the image contrast and 
the temperature contrast would differ unless the object has 
an emissivity value of one. Now, equations to compute 
temperatures at the measurement ROI and the reference ROI 
using the measured ROI intensities will be derived. From 
equations (27) and (31), the temperature at the measurement 
ROI is given by equation (55a), 


o! , 2 
sa Stape Tape — 1 - ©) Tre 
ra| Mawe 


(55a) 


€ 


In equation (55a), all quantities are known or measured 
except the reflection temperature. The reflection temperature 
is derived from equation (35) and may be expressed as in 
equation (55b): 


(55b) 


o { Eft 
2 Wee = 


tape 
reft = 0 . 
Wiape ( l- Efoit ] 
o4 
Trape 


Etape 


The pre-flash temperature at the measurement ROI is given 
by equation (56), 


(56) 


we 425 
fn ster 
won| Wine 


é€ 


In equation (56), all quantities are known or measured 
except the pre-flash reflection temperature, which is given 
by equation (37b). Similarly, the temperature at the refer- 
ence ROI is given by equation (57), 


Ww ri 0.25 (57) 
Frye Toe — (LOT 
W, 
Trop = 


tape 


€ 


Sunilarly, the pre-flash temperature at the reference ROI is 
given by equation (58), 


(58) 


tape 


we 0.25 
ef ot 4 
fare Wyo Tiepe ~ (1 ~ €)Tren 
0 
Tre = 


€ 


Using equation (55a) through equation (58) and equation 
(23) through equation (25), the temperature contrast may be 
computed. Thus, the temperature contrast may be calculated 
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18 
using the ROSE intensity provided by the IR camera. The 
formula assumes that the half-max width at the measurement 
and reference point is at least four times the width of the 
respective ROIs so that the ROI intensity provided by the 
camera is close to the maximum intensity possible for the 
given temperature. 


The accuracy of the estimation of the temperature contrast 
from the measured pixel intensity data is quite sensitive to 
the value of the surface emissivity of the test object. In one 
embodiment, the method assumes that the emissivity of the 
measurement ROI and the reference ROI are the same. In 
another embodiment, the method can either independently 
measure emissivity or use a recommended value. Another 
approach to determine emissivity from the IRFT data taken 
with the foil and the tape on the test object starts with the 
relationship between the pixel intensity and the temperature 
for the foil and the tape. Using equation (36), the method in 
one embodiment can represent the pre-flash pixel intensity 
of the foil and the tape as in equation (59) and equation (60), 
respectively, 


4 4 
Wroit =CoamO (1-€ join) T, eh t€piT, ered ) and (59) 


(60) 


ri a 
On 0 0 
Wrape *C cam Etapel tape +1-Erape) Tren’ )+ 


Equations (59) and (60) may then be arranged as expressed 
in equations (61) and (62), 


wo s,.To! (61) 
tape cys, tape * tape + ne an 
(1 - Exape) (1 — ape) 
4 62 
Wei ee om Efoit Tape 2) 
(leu) "| eon) | 


Equations (61) and (62) may then be combined to express 
equations (63), (64), and (65), 


0 0 
Wiape Wot 


(Sap) (=8pi) 


ape 


4 4 
c ErapeT rape EpoitT eo (63) 
"1 (1 = Eape) (1 = e pot) | 


Weve Wei (64) 


(Sane) (epi) 


Exape 


Cam 
d- Etape) 


Efoil ot 
T, nes and 
d- | saps 


Weoit (65) 


Wiape 2 
(L-f:ape) (1 — Emin) 
E foil ) 


Se Te = 
(a ~Expe) (le pi) 


4. 
Yt 
CcamT ape = € 
tape 


Equation (65) provides a more precise expression for the 
camera constant defined in equation (33). Equations (59) 
and (60) may then be arranged as expressed in equations 
(66) and (67), 
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Wore (1 = Supe JTS : 
2 i CP To eee | oad 
Etape a Etape 
A 
Writ , {a- epiulTon of (67) 
= Cam + Trape + 
Efoil Efoil 
Equations (66) and (67) may then be combined and 
expressed as in equation (68), 
4 4 
Wore Weis eee (L-ftape)Ten (Ll - epi) Tn Ge) 
Sigee sep Sige Efoit : 
Equations (64) and (68) may then be combined and 
expressed as in equations (69) and (70), 
[== : “ é = Exe) (1 a) 2) 
rape Efoit Erape Epil Tea 
0 0 a EY 7 and 
Wrape Wroit ( Frape = Epil Trape 
( — Etape) (dl — foil) a — Erape) a — Epil) 
Were Writ | (70) 
fa ee 4 
Etape foil _ (1 = &tape) (1 - €foit) Tren 
Woe _ Wet Etape Efoit Tore 
dl — tape) a — Eyoit) 


Next, equation (70) may be rearranged and a quantity J may 
expressed as in equation (71), 


Wire — Whit 7) 
a: [= : | 
- Tobe = Wene _ Wit I- (1 = &ape) (1 fs) 
(L-E:ape) (1 — €foit) Etape Efoil 


Further, a quantity B may be denoted as expressed in 
equation (72), 


Wrape Writ (72) 
(a Exp) ( “| 
( Exape Efoil ) 
(L=é:pe) epi) 


Equations (65) and (72) may then be combined and 
expressed as in equation (73), 
B=ChamTape> 


cam* tape 


(73) 


The camera constant is given by the representation in 
equation (74), 


10 


15 


20 


25 


30 


35 


40 


45 


50 


55 


60 


65 


20 
Woe Wha 4) 
(l-Eape) (-efou)} 1 


Cn = 


( Expe Efi 
(l-Epe) (1 -€pou) 


y= 


The pre-flash temperature at the measurement RO] is given 
by equation (75), 


WC say (ET? +H(1-€) Treg). (75) 


Due to some noted observations and assumptions, including 
contact of the tape with the test object, the steady state of 
surface temperature before flash, and the reflection tempera- 
ture not differing more than 5 K from the test object 
temperature before flash, it may be presumed that equation 
(76) holds 


PxT pers 


‘ape (76) 
Equation (76) may then be substituted into equation (75), 
resulting in the expressions of equations (77), (78), (79), and 


(80), 


WC can (Trape H(1-€) Tye )- (77) 
4 a (78) 
re} 
w? = Conte +(l-e) 7 
tape 
W=B(e+1-€)J), and (79) 
we (80) 
e= B 
«o_.. 


Thus, the emissivity of the test object at the room tempera- 
ture can be estimated by an embodiment described herein. If 
the measurement is done at two different ambient tempera- 
tures that are at least 5 K apart, the method may provide 
better confidence in the average emissivity measurement. 
This example illustrates the steps for converting the 
normalized image contrast from the thermography data to 
the normalized temperature contrast using additional mea- 
surements required for determining the reflection tempera- 
ture. In one embodiment, the test object may comprise 
back-drilled elongated flat bottom slots or holes at various 
depths. FIG. 3 shows the IRFT image of a test object that is 
made of reinforced carbon-carbon (RCC) with back drilled 
flat bottom slots and with foil and tape (post-flash time=0.6 
sec). The method in this example is interested in a measure- 
ment on the center slot. The locations of the measurement 
and reference ROIs are indicated in FIG. 3. One piece of the 
wrinkled aluminum foil is placed on the left side of the test 
object and the other piece of an identical foil is placed on the 
right side of the test object. The locations of the foil ROT are 
indicated. A high emissivity adhesive tape is affixed at the 
top of the test object The tape ROI is indicated in FIG. 3, The 
image frame time is extremely close to the peak contrast 
time. The image indicates normalized pixel intensity. The 
image display of the contrast is adjusted by manipulating the 
display rectangle location and size. In this example, the 
measurement ROI is a single pixel (or simply 1x1 pixel). 
The sizes of the ROIs are provided in the Table I. A 1x1 pixel 
in dimensions is equal to 0.023 in.x0.023 in. Normally, ROI 
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cursors and lines have colors for quick identification, but 
colors are not necessary in this description. 


TABLE I 


Pixel Size of Region of Interests 


Region of Interest Pixel Size 
Measurement 1x1 
Reference 3x3 
Foil 2 15 x 15 
Foil 1 15 x 15 
Tape 9x9 


The dashed cross lines of FIG. 3 are used to plot the pixel 
intensity along the lines to locate the peak point and the 
baseline intensity for the selection of the measurement and 
reference ROIs. FIG. 4 shows a plot of pixel intensity along 
the horizontal pixel line taken in a frame at post-flash time 
of 0.6 seconds. Various ROIs, such as the foil, measurement, 
and reference ROI, are also identified in FIG. 4. 

Table II gives the pre-flash measurements and emissivity 
inputs. Note that the representation of Eq. (80) is applied to 
the Table II data (except for the test object emissivity). The 
test object emissivity in this example is estimated to be 0.75 
at the measurement ROI. A significant number of such 
measurements would improve confidence in the value. 
Therefore, the method may instead use an emissivity value 
of 0.78, which is the value published in the literature for the 
object material. 


TABLE II 


Pre-flash Measurements and Emissivity Inputs 


Tye? (K) 300 
tape’ (K) 296 
Erape 0.99 
Gil 0.05 
€ 0.78 
Wrapee (average bit value) 6382 
Wroit’ (average bit value) 6768 
Wie (average bit value) 6480 


FIG. 5 provides plots of the pixel intensity evolutions of 
the two foil ROIs. The vertical dotted line at time equal to 
0.6 seconds indicates the peak contrast time. The foils are 
reflective but specular reflectivity is undesirable. Moreover, 
the test object is also diffusely reflective. The two foils 
provide two independent measurements of the pixel inten- 
sity evolution. The two intensity evolutions agree with each 
other and an average value is chosen. The two foils are on 
either side of the measurement/reference ROI. Thus, an 
average of the two foil intensity evolutions provides an 
estimate of foil intensity evolution in the measurement/ 
reference ROI with better confidence as opposed to choosing 
a single foil intensity evolution. 

FIG. 6 shows the pixel intensity evolution at the tape ROI. 
The flash duration is set to 0.003 seconds. The time between 
the image frames is 0.017 seconds for a frame rate of 60 
frames/second. The intensity decays quite rapidly before the 
capture of the first post-flash frame. The test object reaches 
the maximum temperature at the end of the flash duration, 
but the camera is incapable of capturing the corresponding 
peak intensity due to the slow frame rate. FIG. 6 indicates 
that the temperature rise is higher on the tape ROI than that 
on the foil ROI. The tape ROI is in the area of tape overhang 
and indicates similar intensity evolution as that of the foil 
ROI. 
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FIG. 7 shows the pixel intensity evolution at the mea- 
surement (upper curve) and the reference (lower curve) ROI. 
FIG. 7 indicates the hot spot (higher intensity) of the 
measurement ROI in comparison to the reference ROI. The 
two traces for the measurement ROI and the reference ROI 
eventually merge together as the temperature of the indica- 
tion area and the surrounding area become equal. 

FIG. 8 shows a plot of a simple image contrast, which is 
defined as the difference between the intensities of measure- 
ment and reference ROIs. The simple contrast evolution at 
the measurement ROI, indicates a peak time of approxi- 
mately 0.5 seconds. 

Referring now to FIG. 9, this embodiment begins with 
equation (65) being used to compute the camera constant. 
Next, equations (37a) and (37b) are used to compute the 
reflection temperature evolution. FIG. 9 provides an 
example of a computed reflection temperature evolution. 
The pre-flash reflection temperature is estimated to be about 
305 K, which is 5 K above the pre-flash temperature of the 
test object. Also, the reflection temperature is 2 K higher 
than the pre-flash reflection temperature. When the data for 
this example was taken, the setup was used to acquire five 
other shots within 5 minutes before this represented shot, 
and the time between the last shot and the current one was 
less than a minute. Thus, the flash lamps were warmer than 
the ambient temperature before the flash and caused the 
reflection temperature to be higher than the ambient tem- 
perature. In at least one of the embodiments described 
herein, the method computes the measurement temperature 
evolution using equation (55a). Using equation (57), the 
method also computes the reference temperature evolution. 
The three estimated temperature evolutions are shown 
together in FIG. 10. The upper curve is for the reflection 
temperature. The middle curve is for the measurement ROI 
temperature and the lower curve is for the reference ROI 
temperature. FIG. 11(a) shows the image contrast computed 
using the pixel intensity and one fourth (14) power of pixel 
intensity to illustrate equation (45), which suggests that the 
contrast computed with the fourth power of the temperature 
and the contrast computed with unity power of the tempera- 
ture are approximately the same. FIG. 11(6) demonstrates 
the estimated temperature contrast (upper curve) along with 
the image contrast (lower curve) with e=0.78. Thus, the 
pixel and temperature contrast differ, if the test object’s 
emissivity is less than one. FIG. 11(5) also illustrates the 
difference in the peak contrasts due to lower values of the 
emissivity. FIG. 12 shows an exemplary screen shot of the 
steps of computing the temperature contrast evolution (up- 
per curve) and pixel (image) contrast evolution (lower 
curve) with software implementing the IR Contrast method 
which relates the thermal measurements in terms of the 
normalized anomaly contrast as a function of the frame 
number. FIG. 13 represents the plot of the contrast ratio e'. 
Two plots are shown. One plot (lower curve) uses the 
representation given by equation (52). The other plot is the 
ratio of the two contrasts shown in FIG. 12. The plots differ 
slightly at early instances of time, but the two values agree 
well at later times. The contrast ratio is approximately the 
same as the emissivity of the test object for longer times 
when the rate of change in the surface temperature is 
relatively low. 

In at least one embodiment, the method estimates the 
normalized image contrast based on results from the surface 
temperature simulation of the software sold commercially as 
ThermoCale in order to understand the effect of the reflec- 
tion temperature on the profile of the pixel contrast. In this 
particular example, the method takes a sample simulation 
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temperature evolution at the measurement and reference 
areas. The simulation power input is 1.8E6 W/mA2 with a 
duration setting at 0.003 second. This set-up assumes about 
30% efficiency for converting the electrical power (12 kJ 
source) to the absorbed energy. Depending upon the flash 
power setting, the reflection temperature evolution is 
affected. The method of this example estimates the reflection 
temperature evolution based on the experimental method 
explained in this description. 

FIG. 14 shows the simulated temperature evolutions with 
an assumed reflection temperature evolution. The upper 
curve represents the reflection temperature. The middle 
curve is indicative of the measurement ROI temperature and 
the lower curve is for the reference ROI temperature. The 
normalized temperature contrast is calculated using the 
representation given by equation (51). The method of this 
particular example assumes that the emissivity of the test 
object equals 0.78. Using equation (30), the method of this 
example can compute the image contrast. The results are 
shown in FIG. 15. The upper curve shows the temperature 
contrast. The lower curve shows the image contrast. Using 
equation (27), the pixel intensity can be determined if the 
camera constant is known. Thus, this example illustrates the 
influence of the emissivity change on the pixel contrast. It 
also suggests that changes in the pulse power have a small 
influence on the pixel contrast due to changes in the reflec- 
tion temperature. 

If the method of this particular embodiment assumes that 
one can input the ambient reflection temperature in the 
simulation model, then the afterglow of the heat source 
would he due to differences in the radiation between the 
pre-flash and post-flash reflection temperature. The after- 
glow flux may be represented by equation (81), 


4 
Sposepash=O rep! —Lrep’ )- (81) 


The afterglow flux evolution starts from the first post- 
flash frame. The flash lamp pulse flux evolution is estimated 
for a duration between the start of the flash and the first 
post-flash frame. An example of afterglow flux evolution is 
given in FIG. 16. FIG. 17 shows an example of flash flux 
evolution until the first post-flash frame. The flash lamp 
pulse and the post-flash afterglow evolutions can be added 
to model the compound source flux evolution for the entire 
time of data acquisition. 

FIG. 18 is a flow diagram illustrating methods for con- 
verting an image contrast evolution of an object to a tem- 
perature contrast evolution, for modeling afterglow flux, and 
for assessing the test object emissivity, in accordance with 
some embodiments described herein. Processing begins at 
step 1800 whereupon, at block 1805, a steady state tem- 
perature measurement is made of the tape affixed to the test 
object before the flash or the application of heat. At block 
1810, the transient thermography IR digital video data of the 
object, foil, and tape is collected. The collection of data 
begins before the application of heat. At block 1815, the 
pixel intensity evolutions at the measurement, reference, 
tape, and foil ROIs are measured. At block 1820, the 
normalized pixel intensity ((@mage) contrast evolution is 
computed using equation (20). At block 1825, the camera 
constant is computed using equation (65) and the reflection 
temperature evolution is computed using equation (37). At 
block 1830, the measurement temperature evolution is com- 
puted using equation (55a) and the reference temperature 
evolution is computed using equation (57). At block 1835, 
the normalized temperature contrast evolution is computed 
using equation (23) or equation (51). At block 1840, the 
afterglow flux is modeled using equation (81). At block 
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1845, the test object emissivity is assessed using equation 
(80). Processing subsequently ends at step 1899. 

FIG. 19 is a flow diagram illustrating a method for 
converting a temperature contrast evolution of an object to 
an image contrast evolution, in accordance with some 
embodiments. Processing begins at step 1900 whereupon 
either, at block 1905, the source flux and the reflection 
temperature evolution are assumed, or at block 1910, the 
source flux is based on the heat source setting and modeled 
from the experimental reflection temperature evolution 
using equation (81). At block 1915, the transient thermog- 
raphy surface temperature on a void-like flaw in a selected 
material is simulated. At block 1920, the measurement ROI 
temperature evolution and the reference ROI temperature 
evolution are simulated. At block 1925, the normalized 
temperature contrast evolution is computed using equation 
(23), followed by block 1930 where the normalized pixel 
intensity contrast evolution is computed using equation (30) 
or equation (50). Alternatively, at block 1935, the camera 
constant is calculated from equation (33) and the change in 
pixel intensity evolution at the measurement ROI and the 
reference ROI are computed using equation (28) and equa- 
tion (29), respectively, followed by block 1940 where the 
normalized pixel intensity contrast evolution is computed 
using equation (20). Processing subsequently ends at, step 
1999. 

The normalized image contrast and the normalized tem- 
perature contrast differ for objects with emissivity other than 
one. Therefore, for accurate results, the embodiments and 
examples described herein indicate that the two quantities 
should not be treated as the same. To compare the simulation 
temperature contrast with the measured pixel contrast, the 
method of the embodiments described herein should esti- 
mate the evolution of the reflection temperature and the 
incident heat flux. Ideally, the set of instructions from a 
computer software program that implements the methods 
described herein should model the compound heat source 
flux evolution, which also includes the thermal afterglow. 
The effect of the reflection temperature on, the pixel inten- 
sity should also be accounted for to seek a better estimation 
of the temperature contrast evolution from the pixel intensity 
evolution data. 

As described herein, the reflection temperature evolution 
is established based on data acquired by the IRFT data 
acquisition system. This system acquires readings from a 
test object, a high emissivity tape with known emissivity, 
and a diffuse highly reflective metal foil with known reflec- 
tivity. The method also records the steady state pre-flash 
temperature of the object using a thermocouple (or other 
contact sensor) or an accurate radiometer. The IR datacube 
is recorded using the normal IRFT technique. 

The method comprises four regions of interests (ROIs). 
One region is for the measurement ROI. The second region 
is for the reference ROI. The third region is the foil ROI for 
measurement of the reflection temperature. The fourth 
region is for measurement of the pre-flash temperature at the 
high emissivity tape. 

Using formulas described previously herein, the method 
estimates the reflection temperature evolution. The method 
then computes the temperature contrast from the IRFT data. 
In accordance with at least one embodiment described 
herein, the pixel intensity/temperature contrast ratio has 
been defined in order to relate the temperature contrast to the 
image contrast. 

The embodiments described herein also include a method 
for using the evolution of the reflection temperature to 
model the afterglow flux of the flash source. Using the 
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estimated compound source evolution in simulation soft- 
ware, the temperature contrast evolution may be estimated 
and then the image contrast profiles on the simulated voids 
may be estimated. 

The embodiments described herein include an emissivity 
estimation technique using the same IR camera that is used 
in an IRFT system. This method provides determination of 
the emissivity for the desired thermal wavelength. By using 
the foil-tape technique during the IRFT shot, the transient 
reflection temperature or the reflection temperature evolu- 
tion can be recorded. If an IR camera is programmed with 
the representative formulas for reflection temperature for- 
mulas described herein, the camera can provide the object 
surface temperature directly, even during the IRFT data 
acquisition. The IR camera can also be programmed to 
estimate the object emissivity in real-time by using the 
formulas described herein in combination with the set-up 
described for the foil-tape technique. 

In light of the principles and exemplary embodiments 
described and illustrated herein, it will be recognized that the 
exemplary embodiments can be modified in arrangement 
and detail without departing from such principles. Also, the 
foregoing discussion has focused on particular embodi- 
ments, but other configurations are contemplated, in par- 
ticular, even though expressions such as “in one embodi- 
ment,” “in another embodiment,” or the like are used herein, 
these phrases are meant to generally reference embodiment 
possibilities, and are not intended to limit the invention to 
particular embodiment configurations. As used herein, these 
terms may reference the same or different embodiments that 
are combinable into other embodiments. 

Similarly, although exemplary processes have been 
described with regard to particular operations performed in 
a particular sequence, numerous modifications could be 
applied to those processes to derive numerous alternative 
embodiments of those described herein. For example, alter- 
native embodiments may include processes that use fewer 
than all of the disclosed operations, processes that use 
additional operations, and processes in which the individual 
operations disclosed herein are combined, subdivided, rear- 
ranged, or otherwise altered. 

In view of the wide variety of useful permutations that 
may be readily derived from the example embodiments 
described herein, this detailed description is intended to be 
illustrative only, and should not be taken as limiting the 
scope of the invention. What is claimed as the invention, 
therefore, are all implementations that come within the 
scope of the following claims, and all equivalents to such 
implementations. In the claims, means-plus-function and 
step-plus-function clauses are intended to cover the struc- 
tures or acts described herein as performing the recited 
function and not only structural equivalents, but also equiva- 
lent structures. Thus, while a nail and a screw may not be 
structural equivalents in that a nail employs a cylindrical 
surface to secure wooden parts together, whereas a screw 
employs a helical surface, in the environment of fastening 
wooden parts, a nail and a screw may be equivalent struc- 
tures. 

The invention claimed is: 

1. Amethod for converting an image contrast evolution of 
an object to a temperature contrast evolution, the method 
comprising: 

calculating a measurement region of interest temperature 

change AT; 

calculating a reference region of interest temperature 

change AT; 

calculating a reflection temperature change AT, 
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calculating the image contrast evolution C,,’; and 
converting, using a processor, the image contrast evolu- 
tion C,’ to the temperature contrast evolution C’ 
according to the equation of 


1 
(11st 
C x 1+ TAT SAT) oem 


wherein € is the emissivity of the object. 
2. The method of claim 1, where calculating the image 
contrast evolution comprises: 
measuring a measurement region of interest temperature 
T; 
measuring a reference region of interest temperature T,,. 
measuring a reflection temperature T,g; and 
calculating the image contrast evolution according to the 
equation of: 


ef(1* 1) - (thy - 
(e((T4 - 7?) + (Th 
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wherein € is the emissivity of the object, T° is the measure- 
ment region of interest temperature at time of flash, Toe is 
the reference region of interest temperature at time of flash, 
and T,,,” is the reflection temperature at time of flash. 

3. The method of claim 1, where calculating the image 
contrast evolution comprises: 

measuring a measurement region of interest temperature 

T; 

measuring a reference region of interest temperature T,,,, 

measuring a reflection pixel intensity W,,;; 

calculating a camera constant C'.,,,,; and 
calculating the image contrast according to the equation 

of: 
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wherein € is the emissivity of the object, €,,;, is the emis- 
sivity of a foil, T° is the measurement region of interest 
temperature at time of flash, Toe is the reference region of 
interest temperature at time of flash, and Wroit” is the 
reflection region of interest pixel intensity at time of flash. 

4. The method of claim 3, where the camera constant is 
calculated according to the equation of: 


Woe Wha 
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E foil ] po 


Ce = 
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(l-Eape) (1 - € fit) 


wherein €,,,,. 18 the emissivity of the tape, 1 ae is the tape 
region of interest temperature at time of flash, and Wrape is 
the tape region of interest pixel intensity at time of flash. 
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5. Amethod for assessing an object emissivity, the method 
comprising: 
measuring a pre-flash temperature at a measurement 
region of interest W°; 
calculating a camera constant C',,,,,,; and 
calculating the object emissivity according to the equation 
of: 
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wherein €,,; is the emissivity of the foil, €,,,,. 1s the emis- 
sivity of the tape, Tae is the tape region of interest 
temperature at time of flash, Wri is the reflection region of 
interest pixel intensity at time of flash, and Wishes is the tape 
5 region of interest pixel intensity at time of flash. 

6. The method of claim 5, where the camera constant is 


calculated according to the equation of: 
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7. A method for calculating an afterglow heat flux evo- 


lution, the method comprising: 
measuring a reflection temperature T,,,, and 
calculating the afterglow heat flux evolution S,,,.jasn 


20 according to the equation of: 
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wherein T,,,° the reflection temperature at time of flash and 
o is the Stefan-Boltzmann Constant. 
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